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^ ■ Abstract 

Q^. We propose a three-dimensional electromagnetic current operator within light-front dynamics 

, that satisfies a light-front Ward-Takahashi identity for two-fermion systems. The light-front cur- 
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rent operator is obtained by a quasi-potential reduction of the four-dimensional current operator 
and acts on the light-front valence component of bound or scattering states. A relation between 



> 

o . 

, the light-front valence wave function and the four-dimensional Bethe-Salpeter amplitude both for 



bound or scattering states is also derived, such that the matrix elements of the four-dimensional 
current operator can be fully recovered from the corresponding light-front ones. The light-front 
current operator can be perturbatively calculated through a quasi-potential expansion, and the 
^ \ divergence of the proposed current satisfies a Ward-Takahashi identity at any given order of the 

' expansion. In the quasi-potential expansion the instantaneous terms of the fermion propagator 

are accounted for by the effective interaction and two-body currents. We exemplify our theoretical 
construction in the Yukawa model in the ladder approximation, investigating in detail the cur- 
rent operator at the lowest nontrivial order of the quasi-potential expansion of the Bethe-Salpeter 
equation. The explicit realization of the light-front form of the Ward-Takahashi identity is verified. 
We also show the relevance of instantaneous terms and of the pair contribution to the two-body 
current and the Ward-Takahashi identity. 

PACS numbers: 11.10.St,11.40.-q,13.40.-f,21.45.-v 



1 



I. INTRODUCTION 



A detailed investigation of the electromagnetic (em) properties of hadrons (including also 
nuclei) requests a careful treatment of the current operator for interacting-fermion systems. 
In particular the gauge symmetry plays the well-known essential role, that formally leads 
to the fulfillment of the Ward-Takahashi identity (WTI). Following the seminal paper by 
Gross and Riska [li], in a fully covariant theory, namely in a field theoretical framework in 
its full glory, one should first solve the Bethe-Salpeter (BS) equation and obtain a consistent 



current operator; then through a Mandelstam Q approach, one should evaluate the hadron 
em properties. Unfortunately BS equation can be solved only within approximate schemes 
or for simple examples (see e.g. jsl, 0, O, 0, B])- For instance, an appealing technique 
for solving the BS equation is based on the quasipotential (QP) expansion of the four 
dimensional T matrix (see, e.g. j^), where one introduces an auxiliary four-dimensional 
free Green's function (in jsj it was three-dimensional), hopefully clever enough to allow a 
meaningful truncation of the expansion. Then, one projects the relevant operators onto a 
three-dimensional hyperplane for obtaining a three-dimensional integral equation for both 
scattering and bound states. These three-dimensional states allow to fully reconstruct the 
bur-dimensional BS amplitude solution after applying proper reverse projection operators 
9], |lO|. The method is in principle exact, giving the solution of the four- dimensional BS 
equation if the QP is evaluated without any approximation. As shown in (ll]] for massive 
particles, with a proper choice of the auxiliary free Green's function the QP expansion 
can be put in correspondence with the decomposition of the wave function in terms of Fock 
components on the light-front (LF) hyperplane, i. e. on the three-dimensional space adopted 



in this paper. There are extended reviews 
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14| illustrating the relevant features of 



the LF framework, introduced in a famous paper by Dirac [15j, and the interested reader 
can profit of those. 

The present paper is a part of a program which attempts: i) to solve the four- dimensional 
BS equation for two interacting particles, both bosons and fermions, by using the QP ex- 
pansion and a LF three-dimensional projection onto the LF hypersurface (see o], Q, [3] 
for previous works) and, ii) to develop a procedure for evaluating matrix elements of the 
four-dimensional em current keeping WTI valid, at any order in the truncation. A key role 
is played by the LF reverse projection operator that allows the full reconstruction of the 
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four- dimensional BS amplitudes, without loosing any physics content j^, Q, Q. The price 
to be paid, however, is a rather complicated effective three-dimensional interaction, that 
is derived from the full four-dimensional one. In particular, the effective three-dimensional 
interaction is obtained from the proper projection of the quasi-potential, expanded in pow- 
ers of the four-dimensional interaction present in the BS equation. In conclusion, such 
an approach is useful if the convergence rate of the expansion allows to adopt a truncated 
series, at a nontrivial low order, in the actual calculations. The procedure of LF p^rojection 
and expansion/truncation o 



two-boson case and in Ref. 



jthe three-dimensional interaction is discussed in Ref. ^\ for the 
id] for the two-fermion case. 
By using the two previously introduced ingredients, namely QP approach and LF projec- 
tion, we first investigate the possibility to express the matrix elements of the four- dimensional 
em current operator, for an interacting two-fermion system, through matrix elements of a 
three-dimensional LF current operator between valence states, without introducing any trun- 
cation (see for a formal definition of the valence state). In particular, it will be shown 
that the three-dimensional LF current operator fulfills WTI. This initial analysis will allow 
us to define LF charge operators and the form of the WTI onto the LF hyperplane. After 
introducing a proper truncation procedure of the LF current operator, it will be presented 
a systematical way for constructing workable approximations of the matrix elements of the 
LF current by increasing the order of truncation of the QP expansion. It will be shown 
that the truncated LF current fulfills WTI with the LF charge operators obtained in the 
non truncated case, leaving a complete study of the relativistic covariance under dynamical 
transformations to a future investigation. In order to illustrate our procedure in an actual 
case, it will be analyzed the lowest non trivial truncation of a Yukawa model with chargeless 
boson exchange, in ladder approximation. 

As well-known (see, e.g. Ref. [l^) a peculiar feature of any LF description of fermionic 
systems, is the so-called instantaneous (in LF time) propagation. The treatment of this 
issue makes sharply different the study of the em current for fermionic systems from the 
bosonic case, already analyzed in detail in Ref. [l^, and in what follows this point will be 
often emphasized. 

In general, the construction of a conserved current operator acting on the valence compo- 
nent of the LF wave function is a challenging problem. For two-boson systems Kvinikhidze 
and Blankleider [isl developed gauging techniques for solving this problem. They were able 
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to obtain a current that satisfies a WTI and therefore is conserved. It is remarkable that 
for the two-boson case the current operator derived in Ref. 16| is equivalent to the one 



obtained with gauging techniques 18(]. An earlier perturbative approach to the WTI within 
LF quantization was performed in Ref. fl9|. 

This work is organized as follows. In Sect. II, the four-dimensional expression of the 
em current operator and the associated WTI in the context of BS formalism, that leads to 
current conservation, is briefly recalled. In Sect. Ill, we illustrate the LF time projection 
technique, that allows one to eliminate relative LF time, within a QP approach applied to 
the BS equation for two-fermion systems. In Sect. IV, it is shown the connection of the two- 
fermion valence LF wave function, for bound and scattering states, with the corresponding 
BS amplitude. In Sect. V, the three-dimensional current operator is introduced and the 
fulfilled WTI is discussed. In Sect. VI, we propose a truncated form of the em current 
operator that satisfies the WTI at any given order, and we explicitly show that the current 
conservation is fulfilled by the matrix elements calculated between valence states. In Sect. 
VII, we present the formal construction in an actual case: a Yukawa model with a chargeless 
boson exchange, in ladder approximation. Our conclusions are drawn in Sect. VIII. 



II. FOUR-DIMENSIONAL EM CURRENT OPERATOR 

The general field-theoretical description of an interacting two-fermion system is given by 
the BS equation, where the driving term, the interaction V{K), is assumed to be based on 
the irreducible exchange of bosons. In what follows self-energy corrections will be omitted 
and considered elsewhere. Let us remind that the total four-momentum K is conserved, and 
all BS operators, as V{K), depend parametrically on K (see, e.g. [3]). As well-known, the 
interaction V{K) yields the transition matrix T{K), i.e., 

TiK) = ViK) + VGoiK)TiK), (1) 

and the full Green's function, G{K), i.e., 

G{K) = Go{K) + Go{K)V{K)G{K) = Go{K) + Go{K)T{K)Go{K) , (2) 
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where Gq{K) is the free Green's function of two fermions, with four-momentum operators 
(j=l>2), viz., 

<^o(^) = 77- , (3j 

Kf — mf + le A;| — + 

and k2 = K — ki. For convenience a factor of l/27r has been included in the above definition. 

The BS amphtude dependent upon internal variables, satisfies (the limg^o is under- 
stood) 

G-\K) \^) = 0, (4) 

with appropriate boundary conditions for bound and scattering states 

The em current operator, J''^, corresponding to the interaction V{K) has a free term, J'q, 
and another one, JT/, which depends on the interaction, as dictated by the commutation 
rules between J''^ and the generators of the Poincare group, i.e., 

J'iQ)=Jo'iQ) + JnQ)- (5) 

where JT/ depends parametrically on the 4-momentum transfer Q = Kf — Ki. In particular 



the Lorentz covariance of imposes that (see, e.g. |20|) 



K T{Q') = D'f [WiA-\Kf)] ' J^iQ) D'' [W{A-\K,)] (6) 

where A is a Lorentz transformation, Q' = A"^ Q, Ji{^f) is the total spin of the initial 
(final) system, and is the unitary representation of the Wigner rotation, W{A~^,K), 
corresponding to A. As to the BS amplitude one has 

%mAK, K') = J2 Si'W S,\A)Di,,^,^ [W{A'\ K)] vi/,M;(A:i, K) (7) 



M', 



where k'^ = A~^ ki, K' = A^^ K, and 5*1(2) (A) is the spinorial representation of the trans- 
formation A [ 21[|. 

Notably, the four-current satisfies the following Ward-Takahashi identity (see, e.g., [ll, lZ|) 

Q.J^'iQ) = G-\Kf)e - e G-\Ki), (8) 

where e = ei+'e2 and E,- is the charge operator for the fermion j, with matrix elements given 
by 

{kj\ej \pj) = tjb^ {kj - pj - Q) . (9) 



Since the inverse of the full interacting Green's function is given by 

G-\K)=G,\K)-V{K), (10) 

and 

Q,Jo'{Q) = G,\Kf)e -eG,\K,) (11) 

then the interacting current fulfills the following relation 

Q/.:7r(Q) = e V{K,) - V{Kf)e . (12) 

Once we consider the matrix elements f\J'^{Q) I^E'j), current conservation can be explicitly 
obtained through Eq. (jlj), i.e., 

Q^(vl/^|:r^(g)|vl/,) = (vl/^| [G-\Kj)e -eG-\K,)] |vl/,) = 0. (13) 

III. LF-TIME PROJECTION AND THE QUASI-POTENTIAL APPROACH FOR 
FERMIONS 



Following Ref. lOj, where the LF projection of the BS equation for a fermionic system 
was investigated, here we briefly resume the QP formalism for the LF projection of the 



two-fermion BS equation (see also Refs. 2a, l23|, l2J]), and we introduce a more compact 
and efficient operator notation. 

Crucial for the LF projection is the separation of the fermion propagator in an on-shell 
term and in an instantaneous one, viz. 



k'^-m? + ie k+{k- - k~^ + ^) 2k+ 

where k~^ = {kj_ + m?)/k'^ is the on-minus-shell momentum. The second term in Eq. ([1] 
does not lead to a free propagation in the global time, since the Fourier transform is divergent 
and yields 5{x'^), i.e. an instantaneous (in the LF time) propagation. This term makes the 
treatment of a fermionic system basically different from the treatment of a bosonic one. The 
LF projection is based on the on-shell part of the two-fermion free propagator of Eq. (JSj), 
i.e. 



The role of Go{K) will become more and more clear as the LF integral equations for the 
transition matrix, the valence Green's function, the valence wave function and etc., will be 
derived. It should be pointed out that the instantaneous terms are fully restored into the 
theory through the three-dimensional effective interaction (see below and also Sect. IV). 
Once the full structure of the propagator is taken into account, the matrix elements of the 
three-dimensional em current operator become equal to the corresponding four- dimensional 
quantities, as will be discussed in Sects. IV and V. 

The free Green's functions Gq{K) and Gq{K) are four-dimensional operators which de- 
pend upon the four-momenta of the two fermions. Let us introduce the LF free Green's 
function, go{K), that is a three-dimensional operator which depends upon the LF momenta 
{k^ , ki±) only. It is obtained from the four- dimensional on-shell Green's function by projec- 
tion, i.e. 

go{K) = \GoiK)\ := j dk'fdk^ {k[~\Go{K)\k^) (16) 
= ^e{K-^ -kt)e(kt) 2m,2m,A4fcio.)A^(fc.o„) ^ ^^^^ 

kt{K+ - kt) [k~ - :fcro„, - Kn + ^e) 

where one can choose > without any loss of generality, and A^{kon) = (^^on + /^"^ 
is the positive energy spinor projector. The vertical bar, |, on the right (left) indicates 
that the minus component present in the ket, \k~) (in the bra, {k~\) is integrated out, 
namely th^ou^h the vertical ba. ope.afo. one projects on the LF hypersurface ^ (see 



Refs. [9 



161]). Thus, in this paper the two vertical bars in Eq. f|T6|) do not mean the 
absolute value; they indicate, indeed, the transition from a four- dimensional operator to a 
three-dimensional one. It is worth noting that go{K) yields the free global propagation of 
two on-shell fermions, and its inverse exists in the valence sector, since A_|.(A;io„)A+(A;2on) 
is the representation of the identity in a two-particle spinor space. Furthermore, one 
can express the three-dimensional g^^^K) in terms of the corresponding four- dimensional 
quantity, Gq^{K), as follows (see Appendix A, Eq. (1A6P ) 

g^\K) = %{K) G,\K) n,{K) (18) 

where we define the free reverse LF projection operator (cf the following Eq. (!23l) ) and its 
LF conjugated, as 

Iio{K) := G^{K)\ g,\K) lio{K) := g,\K) \Go{K). (19) 



An explicit expression of Ilo{K) can be found in the Appendix of Ref. lOj. These operators 
have the following properties, as shown by the definitions in Eq. f[T9|) and from Eq. flA6p . 



\MK) = i MK)\ = i 

\Go{K) G^\K) Ilo{K) = I no(ir) G^\K) G,{K)\ = I (20) 

where I is the identity in the LP three-dimensional space. 
The LP free state, |0o)) is a solution of 

g^\K) |0o) = (21) 

then from Eq. (fTSi) one has 

MK) G,\K) no(ir) |0o) = (22) 

The last equation leads to the following transformation property between the non-interacting 
BS amplitude and the corresponding LP three-dimensional state, viz 

|vl>o) = MK) |0o) (23) 

Indeed, one obtains the eigenequation for |\I^o) by applying Gq ^(f^) to Eq. fl23|) and per- 
forming the understood lim^^o- Then, on the rhs of Eq. (!23|) one can use Eqs. (fT9l) and 
( ]A3[) (that holds for any e) and can shift lime_»o to g^^{K) |0o); namely one has 

Go'(i^) |^o) = (24) 

Purthermore, one gets two expressions of the inverse of Eq. (l23l) by using Eqs. ( l20l) and 
(El, viz 



ll^o) = |0o) (25) 

and 

\Go{K) G^\K) = |0o) (26) 



In Eq. ( I26ll the understood lim^^o affects the whole Ihs. Notably, from the uniqueness 
of the solutions of Eq. fl21|) and applying Appendix A of Ref. [lo|, one can demonstrate 
that for any given solution |\I'o) of Eq. fl24l) one has a LP wave function obtained by Eq. 
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(126|) . that in turn fulfills the eigenequation (l2Ti) and yields the initial |^o) through Eq. (l23l) . 
Namely, there is a one-to-one relation between a given four-dimensional BS amplitude l^'o) 
and the corresponding three-dimensional |0o)- The non-interacting operators Ilo{K) and 
Ilo{K) connect three- and four-dimensional quantities. In the next Section the corresponding 
interacting operators will be given. 

The QP formalism makes use of the four-dimensional auxiliary Green's function Gq{K) 
defined by 

Go{K) := Go{K)\ g,\K) \Go{K) = Uo{K)go{K)Uo{K) = 

= Go{K)\ Uo{K) = Uo{K) \Go{K) , (27) 

This operator has the following useful properties 

\GoiK) = \GoiK) Go(K)| =Go(K)| . (28) 



Let us apply the QP formalism to the four- dimensional transition matrix (see Refs 



10|), i.e. 

T{K) = W{K) + W{K)Go{K)T{K) = W{K) + T{K)Go{K)W{K), (29) 
where the effective interaction W{K), in turn, is a solution of 

W{K) = V{K) + V{K)Ao{K)W{K) , (30) 

with 

Ao{K) := Go{K) - Go{K) . (31) 

The four-dimensional quantity Ao{K) clearly contains the instantaneous terms and has the 
following properties 

71+72+ Ao(K)|=0 7i+ 72+ |Ao(K) = (32) 

It is worth noting that in Eq. ( 129|) the driving term V{K) is substituted by W{K), in order 
to increase the rate of convergence of the iterative solution of the transition matrix (for 
numerical studies of the convergence see, e.g., [sl, Q]). In a similar manner, Eq. (!30l) can be 
solved by iteration, obtaining 

oo 

W{K) = Y,W.{K), (33) 

i=l 



with Wi{K) = V{K) [Ao{K)V{K)f-\ The physical meaning of the convergence of this 
series will be discussed in Sect. VI. 

The three-dimensional LF transition matrix can be introduced through (see also [lo| ) 

t{K) = Uo{K) T{K) Uo{K) (34) 

Then, by using Eq. (I29l) . one has 



tiK) = noiK) \w{K) + W{K)GoiK)TiK)\ UoiK) = 

= w{K) + w{K)go{K)t{K) = w{K) + w{K)g{K)w{K) , (35) 

where the three-dimensional driving term, w{K) is obtained from the four- dimensional in- 
teraction W{K) according to 

w{K) := Uo{K)W{K)Uo{K) . (36) 

and the three-dimensional interacting Green's function is given by 

giK) = go{K) + g,{K)t{K)g,{K). (37) 

Note that i) the positions of Hq^K) and Iio{K) in Eq. ( IMll . leading to the integrations over 
the external minus components, straightforwardly indicate that t{k) is a three-dimensional 
quantity (for comparison see Eq. fl27|) . where the four-dimensional quantity Gq{K) is defined, 
and the integrations are on the internal minus components); ii) the effective interaction w{K) 
contains the coupling of the valence sector to the higher Fock-state components of the wave 



function through W{K) 11 | 



Finally, the four-dimensional T{K), can be rewritten as follows 

T{K) = W{K) + W{K) Uo{K) g{K) Tlo{K) W{K) (38) 

Such an expression is useful to discuss the relation between the four-dimensional BS ampli- 
tude for a bound state and the corresponding three-dimensional valence wave function (cf 
Sect. IV). 

The LF interacting Green's function, or resolvent, can also be written as a function of 
w{K) as follows: 

g{K) = go{K) + g^{K)w{K)g{K) = go{K) + g{K)w{K)g^{K). (39) 

10 



It should be pointed out that g{K) is the Fourier transform of the global propagator of 
two interacting fermions between two LF hypersurfaces with given a;"'"'s. Its inverse can 
be easily related to the four dimensional G~~^{K) by using Eqs. (IT5]1 . (15^ and (|5U]) . As a 

matter of fact, one has 

g-\K) = g,\K) - w{K) = UoiK) [G,\K) - WiK)] UoiK) 
= Uo{K) [G,\K) - V{K) (1 + Ao{K) W{K))] Uo{K) 

= lioiK)G-\K) [1 + Ao{K) W{K)] Uo{K) - %iK)G^\K)Ao{K)W{K)UoiK)iAQ) 



The last term is vanishing, as can be shown by using Eqs. f[T9|) . flA4p and the second 
relation in Eq. (!32|) . Finally one obtains 

g-\K)=n,iK)G-\K)U{K) (41) 

where 

U{K) = [1 + Ao{K) W{K) ] Uo{K) (42) 

is the interacting LF reverse projection operator. From analogous steps one can obtain the 
corresponding LF conjugated operator, given by 

U{K) = Uo{K) [1 + W{K) Ao{K) ] (43) 



It will be very useful in what follows to note that by using Eqs. flA3l) . flA4l) . fl20l) and fl32|) 
one has in the three-dimensional space 

\Go{K) G^\K) n{K) = I T[{K) G^\K) G^{K)\ = I (44) 

The solution of the following three-dimensional equation, 

g-\K) 10) = [g,\K) - w{K)\ |0) = (45) 

with appropriate boundary conditions for bound and scattering states, is the valence compo- 
nent of the LF wave function jisl. It turns out that the full complexity of the Fock-space 
comes through the effective interaction. However, the truncation of the quasi-potential (see 
Eq. ( l33l) ) limits the number of Fock-components involved in the construction of effective 



interaction to be used to obtain the valence wave function 
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ll|. Then, one could argue that 



the convergence rate of the QP expansion is related to the smallness of the probabihty for 
the higher Fock-components. 

Inserting Eq. (1411) in the eigenequation (I45p one has 

g-\K) 10) = no(ir) G-\K) Ii{K) = (46) 

that leads to the following relation between the three-dimensional valence component and 
the four- dimensional BS amplitude 

|vl/) = n{K) 10) (47) 

Note that |\E') fulfills Eq. (|4]) (as can be seen by using the results in Appendix A of flo|). 
Furthermore, by applying Eq. (jH]), one gets 

\G,{K) G-,\K) \^) = 10) (48) 

In the next Section more details will be given about the one-to-one relation between the 
four-dimensional |\1/) and the three-dimensional |0). 

Finally, let us remind that, the on-mass-shell matrix elements of T{K), which define the 
two-fermion scattering amplitude are identical to the ones obtained from t{K) 



IV. THE LF VALENCE WAVE FUNCTION AND THE BS AMPLITUDE 

In this Section, we will analyze the interacting operator Il{K), Eq. P2|) . that generates 
the full BS amplitude, for bound and scattering states of two-fermion systems, starting from 
the corresponding valence wave functions. It should be emphasized the key role played 
by Il{K), when the matrix elements of a four-dimensional operator acting on the BS am- 
plitudes are considered. In particular, the LF reverse projection allows to express those 
matrix elements in terms of matrix elements of effective operators acting on the valence 
wave functions. 

The relation between the BS amplitude for a two-fermion bound system, I^&b), and the 
corresponding valence wave function has been derived in jioj], and reads as 

\^b) = Go{Kb)W{Kb) IIo{Kb) \(Pb) ■ (49) 
This can be obtained through the analysis of the poles of Eq. (|38|) . 
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Since the valence wave function is the solution of 



\<Pb) = g^{KB)w{KB) , (50) 

the vanishing quantity Gq{Kb)\ {gQ^{KB) — w{Kb)) = can be added to Eq. fHOl) in 
order to recover the LF reverse projection as given in Eq. (H2l) . i.e. 

\^b) = [1 + Ao{Kb)W{Kb)]Uo{Kb) \<I)b) = 11{Kb) \<Pb) ■ (51) 

Therefore any {"^b) can be generated through II{Kb) from the corresponding valence wave 
function \4>b), eigensolution of Eq. fHSjl . Moreover, from Eq. fj48|) one has 

\Go{Kb) G^\Kb) \^b) = (52) 

Let us now show that the operator n(i^) connects BS amplitudes and valence wave func- 
tions for scattering states, as well. 

The BS amplitude for scattering states satisfies the four- dimensional Lippman-Schwinger 
type inhomogeneous equation: 

|^+) = \^,) + G,{K)T{K)\^!o) = [l + G,{K)T{K)]Uo\ct>o) , (53) 

where we have made use of Eq. fl23l) . that univocally relates |\&o) to |<^o)- In the three- 
dimensional space the scattering solution of Eq. (l45l) is given by 

\^^) = [l+giK)w{K)]\<Po) 

= [go{K)+g{K)w{K)g,{K)]g,\K)\ct>,) = g{K)g,\K)\<P,) , (54) 

which implies the formal identity gQ{K)g~^{K)\(l)'^) = \(f)o). Then, by using Eq. fIBip . one 

gets 

|vl/+) = [1 + Go{K)T{K)] MK)\^o) = 

[l+Go{K)T{K)]Uo{K) goiK) g~\K) = 

= n(ir)|0+). (55) 

From Eq. (I48l) . the valence component of the LF wave function can be obtained directly 
from the BS amplitude, i.e. 

\Go{K)G,\K) |vl/+) = 10+) (56) 
13 



It is worth noting that the operator \Go{K)Gq (K) cuts the instantaneous terms and 
projects onto the LF hyperplane (through the k~ integration), while Il{K) reconstructs 
the full structure of the four- dimensional BS amplitude through the instantaneous terms 
contained in Go (see Aq) and W. Moreover, the instantaneous terms affect the effective 
interaction w that determines and (cf Eq. fH5l)). 

We should point out that the normalization of the BS amplitude for a bound state {27!] ex- 
pressed in an operatorial form, Af{K), can be mapped onto an expectation value of a proper 
three-dimensional operator viz., (\E'|A/'(-ft')|^) = {(t)\U{K)Af{K)Il{K)\(f)) . Note that on one 
side, the BS amplitude normalization corresponds to the sum over the probabilities of each 
Fock component in the full LF wave function [17], and on the other side the full complexity 
of the Fock-space is summarized in the three-dimensional operator Il{K)^/{K)Il{K). 



V. LF WARD-TAKAHASHI IDENTITY 



Once we have the relation between the 4-dimensional BS amplitude and the 3-dimensional 
LF valence wave function (see Eq. (I47l) ). the LF em current operator can be obtained from 
the matrix element of the four- dimensional current. As a matter of fact, one has for both 
scattering and bound states 

{^jU'^iQ)]^,) = {<PfmKf,K,m , (57) 

where the three-dimensional LF current operator, acting on the valence wave functions, is 
defined as follows 

f{Kf, Ki) := T[{Kf)J^^{Q)Ii{K,) . (58) 

Using Eqs. (jUj) and fH5]l . one can put in evidence the dependence of the LF current operator 
upon the effective interaction W{K), viz 

/(ir^, K;) := T[o{Kj) [1 + W{Kf)^^{Kf)] J^^^Q) [1 + /\o{K,)W{K,)] Il^{K,) . (59) 

It is worth noting that Eq. flS^ generalizes to scattering states the expression for the three- 
dimensional current already derived in Ref. 1^ for bound states. It should be pointed 
out that is Poincare covariant (cf Eq. ([6]) for the Lorentz covariance). Since all the 
matrix elements of the Ihs are properly related through the Lorentz transformations, the 
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investigation of the covariance properties of the operator j^, within the full theory, does 
not represent a stringent question in view of the equality in Eq. fl57|) . But, given the 



id- 



following development of the truncated approach, where the full covariance is broken, it is 
important to investigate the transformation properties of j^. As discussed elsewhere j25|, we 
can anticipate that the covariance with respect to the kinematical transformations can be 
demonstrated after introducing new factors in the vertical bar operation (cf also the factor 
h in Ref. 111). 

In what follows, our task is to find a suitable definition of the LF charge operator that 
allows to write the four-divergence of the LF current operator in terms of the inverse of the 
Green's functions g^^{Kf) and g^^{Ki). Such an investigation is fundamental to obtain 
the WTI for the truncated LF current (see the next Section). It should be pointed out 
that a similar analysis has been already performed for bosonic systems in 

In order to find the LF charge operator, the four-dimensional divergence of the LF current 
can be written as follows by using Eqs. ( !58l) and ^ 

Q^f{Kf,Ki)=T[{Kj) [G-\Kf)e-eG-\Ki)\ n(i^,) , (60) 

with e = ei + 62- 

By applying Eqs. flB3|) and flB4l) of Appendix B, one gets 

Q^f{Kf,Ki) = g-\Kf) \G^{Kf) G^\Kj) e Ii{Ki) -T[{K j) e G^\K,) Go(ir.)l 9~\K,) = 
= g-\Kf) - Qfp g-\K,) (61) 

where left and right LF charge operators have been introduced. Such operators are defined 
as follows 

Qi^ = \Go{Kj) G^\Kj) eUiK,) = \GoiKj) G^\Kj) eUom (62) 

Q^p = n{Kf) eG^\K,) Goim = Uo{Kf) e G^\K,) Go{K,)\ (63) 

where Eqs. ( IA3I) . flA4p and (!32l) have been used, as well as the absence of a Dirac structure in 
the operator fT. It is very important to note that the LF charge operators are non interacting 
operators acting on the two-fermion spinor space. From Eqs. (lA4p and (lASp . we can obtain 
the explicit expression for the left LF charge operator for particle 1 

Q^LF = A+(klor,)^lt eiLFA+(A?ion)A+(A?2on) • (64) 

15 



where the notation 'enF indicates the three-dimensional LF counterpart of the operator e, 
Eq. ([9]), with matrix elements given by 

k[^\eiLF\kt, ■■= eiS {k'+ - kf - (k[^ - h± - Q±) , (65) 

The corresponding right operator is 

Q^LF = Kiklon) eiLF^lt^+{klon)A+{k2on) ■ (66) 

k{ 

The operator '-)^m/k~^ when sandwiched between LF spinors gives the normalization con- 
dition. Let us remind that such an operator is the one particle free charge operator of our 
model, since we do not include fermion self-energy. This also indicates the ingredients that 
ought to be considered when the full problem with self-energy insertions is aimed. We will 
not discuss further this issue here, which is left for a future study. 

From Eq. (l6T]) . current conservation straightforwardly follows by taking the matrix ele- 
ments between three-dimensional interacting states that are solutions of the wave equation 
(ll5ll and noting that the left and right charge operators do not contain any le dependence. 

By multiplying both the left and right hand sides of Eq. ( 1611) by g{Kf) and g{Ki), 
respectively, one gets 

Q,giKj)riKj,KMKi) = Q^MK^d ~ 9iKf)Q^^ , (67) 

which corresponds to the LF projection of the 5-point function with instantaneous terms 
cut from the external fermion legs. 

In conclusion the three-dimensional LF current given by Eq. fl58|) . acts on the LF valence 
wave functions and fulfills WTI (see Eq. (EZl))- 

VI. WARD-TAKAHASHI IDENTITY FOR THE TRUNCATED LF CURRENT 

In this section, to obtain a workable approximation for the LF current operator that still 
fulfills WTI, we will analyze the consequences of a proper truncation of the QP expansion in 
the formal solution (Eq. fl55]) ) of Eq. fl5Ul) . As already anticipated in the previous Section, 
the LF charges and the general form for the LF WTI (cf Eq. (I6T1) ). will represent our 
fundamental ingredients. To simplify notations, we will not show explicitly the parametric 
dependence on total momentum K in the operators. Therefore, it is understood that the 
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operators on the left of the charge operator depend upon Kf, while the ones on the right 
depend upon Ki. 

A naive substitution of W in the current operator (15^ . by the truncated quasi-potential 
expansion 

n 

iyW = ^iy. with Wi = V[AoVy-^ = [VAoY-^V (68) 

i=l 

namely by the truncated iterative solution of Eq. (1501) . does not lead to a conserved three- 
dimensional current. The same problem has been met and solved for two-boson systems 
in Ref. In this Section we will follow the same procedure, but applied to the case of 
two-fermion systems. In the next Section an actual application for the Yukawa model will 
be illustrated. 

Let us first consider the conserved current, for the case n = 0, i.e. the current without 
interaction. In this case one has 

:= HoJo'^no = g^'\GoJ,^Go\9^' . (69) 
From Eqs. ([69]), ([II]), ([62]) and ([63]), WTI for n = reads: 

Q%^'^ = 9o' I Go Go ' e Ho - Ho e G^' Go\ g,^ = 

= ^7o'^QiF-QW- (70) 

The matrix elements of Eq. ( jTOl) should be taken between solutions of gQ^\4>o) = 0. 

The matrix elements of between free particle states, as obtained explicitly in Ap- 
pendix C together with the corresponding WTI, are given by 

(fc'+fc'i^ kthi.) = -te {kt) (A;'+) e {Kt - kt) e {k+ - k'+) 

A+(A;U)7rA+(A;i,„)A+((K, - k,U) ^\~ {k'+, k[^\e,,LF\kth±) + 1^2 (71) 

From Eq. (I59l) and cutting at the first order the effective interaction (note that W^^^ = V) 
one has 

where the first order contribution is given by 

/(Dm = Ho [J^ + VAoJo^ + Jo^AoV] Hq = 

= j^(o)'^ + Ho [J/ + VAoJ^' + Jo^AoV] Ho , (73) 
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Only j'^(^)'^ is a conserved current operator in the corresponding valence sector. Indeed, as 
shown in details in Appendix D, one obtains the following WTI 



QX^'^ = 91-' QIf - QIf 91-' 



(74) 



where gi~^ = go~^ — w*-^^ and w^-^^ = Hq V Hq. 

The matrix elements of Eq. flTll) should be taken between solutions of (7jf^|</)i) = 0. In 
Sect. VII, we will discuss in detail the first-order LF current for the Yukawa model in ladder 
approximation, and we will show that it is essential to include the instantaneous terms, 
coming from Aq to obtain WTI for n = 1. 

For the general case, where the interaction appears up to n > 1 times in the LF current 
operator, we write 



J 



c{n)fj, 



n-1 



i=l 



Tin . 



n-1 j-1 

+WiAoJ^ + J'^AoWi) + W^AoJ^AoWj 

j=2 i=l 
n n—1 

WiAoXAoWn-^ + J2 W.AoJi'AoWn-l-i 

. i=0 i=Q 



(75) 



where it has been formally defined VTo^o = AqWo = 1- It is worth noting that Eq. (!75|) 
contains power of the interaction up to the ra*'' order, since J'q is 0{V^) and ^7/ is 0{V^). 

As demonstrated by induction in Appendix E, this truncated current operator satisfies a 
WTI given by 



/nM.-c(n) _ -1 qL _ qR -1 
^ Jn — 9n ^LF ^LF 9n 



(76) 



where ^ = 9o ^ — w^"^ and w^'^'^ is the truncated effective interaction, viz 

n 

'("^ = UoW^^^Uo = Y ^oW.Uo. (77) 



i=l 



The matrix elements of Eq. (17^ should be taken between solutions of g:^^\(j)n) = 0. 

Thus, we conclude that the LF electromagnetic current operator j^''"^ is conserved at 
any given order n of the QP expansion. In the limit n —>■ oo the truncated conserved 
current becomes the full current operator of the model, if the QP expansion converges (for 
numerical studies of the convergence see e.g. sl, [3]). Finally, the effects of the truncation in 
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the QP expansion, and correspondingly in the Fock-space expansion, on the transformation 



25|. 



properties under the Poincare group will be investigated elsewhere 

It should be pointed out that the last line of Eq. (I7SI) suggests a possible diagrammatic 
representation of the current (cf the next Section for the first-order case in a simple example), 
through the observation that the contribution to the quasipotential at a p*'^ order, Wp, 
can be decomposed in Wp = Wi Aq Wp-i ioi p > i > (see Eq. (!68l) ). Indeed one 
can construct a picture of the contribution between square brackets in the last line of Eq. 
(!75|) . by starting with all the possible decompositions for Wn and Wn-i, then breaking the 
subtracted propagation represented by Aq and inserting the current j7o and J}, respectively, 
i.e. Ao AoJo(/)Ao- 

As shown in detail in the example of the next Section, the subtracted Green's function 
Aq eliminates the presence of disconnected diagrams, due to the global, free-propagation 
of the two-fermions between the photon absorption and the exchange of a boson from the 
interaction. At the same time, it allows for the photon being absorbed i) instantaneously (in 
LP time!), or ii) by higher Pock-state components, including a particle-antiparticle pair. In 
particular, one can easily see that the disconnected diagrams generated from the gauging on 
the external legs, i.e. diagrams that are already taken into account through the interacting 
LP wave function, are eliminated by the presence of the subtraction in Aq, avoiding double 
counting. 



VII. LF CURRENT AND WARD-TAKAHASHI IDENTITY IN THE LADDER 
YUKAWA MODEL AT THE FIRST ORDER 

In this Section we explicitly derive the current operator f^^^^", see Eq. (1731) . for a two- 
fermion system interacting through chargeless bosons, within the Yukawa model in ladder 
approximation. This will illustrate in an actual case how our procedure works to obtain the 
conserved current operator at the first nontrivial order, namely n = 1. The Yukawa model 
is defined by the following interacting Lagrangian (with no derivative coupling): 

Ci = gij^V^ipiaa + 5'V'2r2 V'20-Q, (78) 

where g is the coupling constant (that should not be confused with the LP Green's function 
g{k)\), ipi and tp2 are fermionic fields, the corresponding particles have masses mi and m2, 
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and charges ei and 62, respectively. The field a a corresponds to a chargeless boson with 
mass yU. The index a represents Lorentz components in the case of a vector or tensor field. 



where the symbol ® takes distinct the vertexes corresponding to the fermions 1 and 2. 

We choose this simple example to show that an interacting LF two-body current can be 
generated even starting from a four- dimensional free current, Jq . It should be pointed out 
that in ladder BS approximation the consistent four-current is the free one, i.e. Jq , 

and the four- dimensional WTI is fulfilled because the commutator [e^, V] trivially vanishes. 
Such a cancellation is due to the fact that the exchanged momentum in the potential just 
depends on the spectator particle momentum, that remains unchanged in the transition 
from the initial state to the final one (cf. Ref. \\ and Fig. [1]). Since we need all the four 



FIG. 1: Diagrammatic representation of the matrix elements of the 4-dimensional operator 
GoJqGqVGq, (cf. Eqs. (|Fip and ()F2p ). Dashed line: four-dimensional interaction V. The 
full dot represents the free current. 

components of j'^(^)^ to construct the WTI, Eq. (17^ . in what follows we explicitly evaluate 
all the components of the truncated em current, by projecting the four- dimensional current 
onto the LF hypersurface (i.e. by integrating on k~). 
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Let us start with the matrix elements of the free current operator, given by 

{ki\J^mVi)= -2ne,j^6\h-p,-Q) [(^ f - - m2)] 

+ [l^2,h^Kf~h,p,-^K,-p,] , (80) 

where = — . The factor (— 27r) is introduced in the current operator to make it 
compatible with the free Green's function, see Eq. ([3]). 

The first-order current operator for the interacting two-fermion system in this example 
is given by (cf Eq. fl73|) ). 
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c(i)M = _^-c(o)M ^ no [V^,J^^ + J^^,V] Ho , (81) 



where the zero order LP current operator is the LP free current, see Eqs. f l69|l and fl7T|) . 

The contribution of the interaction to the LP current operator in lowest order comes 
from two-body irreducible amplitudes given by the second term in the r.h.s of Eq. (1731) with 
= Jq , since we are adopting the ladder approximation. Using Eq. fl3T]) for Aq we write 
that: 

^.c(i)^ - = Ho y Go Jo^ Ho - w^^^g,3'^'^^ + Ho J,^ G.VIl,- f^'^''goW^'\ (82) 

where w^^^ = Uq V Uq is the three-dimensional effective interaction, with matrix elements 
for a total momentum K given by (see Appendix C in [lo| ) 

{k['-k[^\w^^\K)\kiki^) = % {%gf e{k'+)e{ki)e{K+ - k[+)e{K+ - k+) 



X 



[K^ - kt) ^ e{kt- k'^) 



{k^ kf) {^K k^^,^ A;2o„ k„on + 'i^) {kf k'^) {^K k[^^ /c2o„ + ka-on + 'i^) 
X A+(A;;„jr?A+(A;io„)A+(A;;,jr2aA+(fc2on) • (83) 

where k2 = K — ki, 

, _ _ fcix + ml _ k[\ + mj 

_{K±^k^^l±tl 

and the quantities corresponding to fermion 2 easily follows. As shown in Appendix P, both 
the first term and the third one in the r.h.s. of Eq. (182!) contain a LP two-body reducible 
part (already taken into account when is apphed to the valence wave function), that 
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is canceled by the second and the fourth ones, respectively. In particular (see Ref. |10|). 
the LF two-body reducible contributions are generated by the effective interaction used to 
obtain the valence wave function at the first order in the quasi-potential expansion of the 
ladder BS equation. The essential role played by the reducible terms was stressed in 



in a calculation of higher Fock states contributions to the Generalized Parton Distribution 
of pion. 



16j, comes from the instan- 



The difference with the boson case, that was recently studied 
taneous terms present in the quasi-potential expansion of the fermionic current operator, as 
we will show in detail. 

In the following we report the matrix elements of the first-order current operator relevant 
for the em processes in the spacelike region, as obtained in Appendix F. 



A. The term UoX{Q)GoVUo 

In Fig. 2, it is diagrammatically illustrated the third term in the r.h.s. of Eq. (I82l) . 
go^\GoJ^Q{Q)GQVGQ\ Qq^. It is worthwhile to note that this term contains the relevant 
pair production contribution, for > (cf. diagram (b) in Fig. 2). Differently, in the first 
term, gQ^\GoVGoJ'Q{Q)Go\ gQ^, there is no contribution from the pair production by the 
virtual photon, since the conservation of the plus momentum component always requires a 
positive plus component of the intermediate fermion momentum and then the initial state 
cannot contain a photon and a pair (cf. Fig. 3). 

Performing analytical integrations over A;f and k[ by Cauchy's theorem with the condi- 
tions > and > (see Appendix F), we get the six contributions that appear in 
Fig. O We observe that the limited number of LF time ordered diagrams are essentially a 
consequence of the trivial vacuum structure, due to the conservation of the total positive 
plus momentum. Now, let us discuss in detail the irreducible diagrams of Fig. O 

Diagram (a) gives a two-body current with intermediate three-body virtual state propa- 
gations before and after the photon absorption, in the kinematical range pf = k'^ — > 
and kf > {k'^ — Q^). Diagram (b) represents the contribution to the current of the vir- 
tual photon decay in a pair. The kinematical condition in this case is Pi = k'^ — < 
(note that —Q^ < Pi < 0). Diagram (c) represents a contribution from the instanta- 
neous term of the fermion propagator, which is non vanishing in two kinematical regions: i) 
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FIG. 2: LF time-ordered diagrams representing the matrix element of the current operator obtained 
from IIqJ'q{Q)GoVIIo, i.e. the third term in the r.h.s. of Eq. ([82]) . The reducible diagrams, (d) 
and (e), are canceled by j^^^^^^gow^^\ 

pt = k[+ - Q+ > and kf > {k'+ - Q+), ii) p+ = k[+ - Q+ < 0. Diagrams (d) and (e) 
are two-body reducible terms which are canceled by j^^^^f^gow^^^ in Eq. ( |82l) . Diagram (f) 
represents a contribution of the instantaneous term for p'f = k'^ — > k^ > 0. 

Defining 6 = 0{kf) 9{k^')6{Kf — kf) 0{K^ — kf), the diagram (a) is given by (see Eq. 

m) 

e {kt - pt) {pt) 



{kt-pt)pt A^^'^^-Af")" 



xA+(A;U)7f A+(pi„„)r^A+(fci„OA+((ir; - k[),^)T2^A+{{K, - k,)on) ■ (85) 
The four-vectors A^^"^ and A^^^'*'' are the following combinations of the four-momenta 

A^^"'' = Kf - k[^^ - {Ki - ki)on - {ki - pi)on + ie 

Af = Ki- Plan - {Ki - ki)on - {h - Pl)on + , (86) 

The minus components yield the three-body (2 fermions and the exchanged boson) global 
propagation in the final and initial states, respectively. 

The matrix element of the current operator (b) (virtual photon decay in a pair) is (see 
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Eq. dFTH)) 

X A+(fcU)7!'A+(p,<,„) r?A+(fci„)A+((A', - *;)„)r2„A+((A', - ki)„,) , (87) 
where the combination A^^*"* of four momenta is 

) =Q- k[,^ + Plan + ^£ . (88) 

The instantaneous term (c) (see Eqs. ( 1F9P and ( 1F15I) in Appendix F) can be expressed 



as 



1 {kt - pi) 



X A+(A;;,j7f7+r?A+(A;io„)A+((i^/ - A;;),„)r2aA+((K, - A;i),„) (89) 

Finally the contribution from the instantaneous term for k'i' — > kf, represented by 
diagram (f), is (see Eq. (]F11|) ): 



2 ' ' Pt (j9+-A;+)Af^)- 

X A+{k[^Mltm+ihon)A+{{Kf - k[U)r2'^A+m - h)on) , (90) 

where the denominator is the minus component of the following four-momentum 

Af ^) = K, - fci,„ - {Kf - k[),^ - (pi - A;i)„„ (91) 

and yields the intermediate propagation of the three-body system composed by two fermions 
and a boson. 

B. The term TioVGoJ^{Q)Ilo 

Let us now discuss the first contribution in Eq. (!82|) . go^\GoVGoJ'Q{Q)Go\gQ^, depicted 
by the diagrams (g) to (i) in Fig. [3l It is worth noting that the pair contribution is not 
present due to the conservation of the plus momentum component. For the sake of simplicity, 
the two reducible diagrams, which are canceled by in an analogous 

way as we have already discussed for diagrams (e) and (d) in Fig. [21 are not shown in 
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FIG. 3: LF time-ordered diagrams representing the matrix element of the current operator obtained 
from IIqVGoJ'q{Q)IIo, i.e. the first term in the r.h.s. of Eq. ([82]) . The reducible diagrams, canceled 
by w^-^^ j'^^'^^^go, are not shown. 

Fig. ([3]). As discussed in detail in Appendix F, by analytical integrations over ki and k[ by 
Cauchy's theorem one can express diagram (g) as (see Eq. (]F21|) ) 



X A+(*;„.)r?A+(p;„,)7f A+(*:i„„)A+((/0 - k',)^)r2cA4{Ki - h)„„) , (92) 
where the combinations A^^^^ and A^^^^'' of four momenta are 

Af ^) = Ki~ k,on - {Kf - k[)on - {k[ - p[)on + ^e = Af ^) . (93) 

since k[ — p[ = pi — ki. 
The instantaneous diagram (h) is given by (see Eq. ( 1F23I) ) 



2 (fe'i^-p^K-^Af^)- 
xA+(A;U)F?7i+7rA+(fcio„)A+((i^/ - k[),^)V,^k+{{K, - k,U) , (94) 

while the instantaneous contribution in the region p'^ > k'^ > 0, illustrated by diagram (i) 
in Fig. da]), is (see Eq. flFM]) ) 

^ [Pi - K)vi 

xA+(fcU)F?7+7fA+(A;io„)A+((iry - k[)on)T2o.K{{K, - ki)on) , (95) 
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where the three-body denominator is the minus component of 

Closing the discussion on the first-order current, we refer the reader to Appendix G for an 
exphcit check of the current conservation. The detailed derivation of current conservation 
allows a deeper understanding of the explicit effects of the formal manipulations used to 
obtain a conserved LF operator, consistent with the effective interaction at each given order. 
This calculation also shows the essential role played by the instantaneous terms in the 
cancellation of terms from the two-body current (a) and (g), that otherwise break current 
conservation. In this respect it is useful to consider that current (i), as current (c), yields 
contribution in both the kinematical regions i) = k'^ — Q^ > and ii) > k'^^ — Q^ = , 
as easily obtained from the kinematical constraint p'^ = kf + > k'^ leading to ki > 

The peculiarity of the instantaneous terms in the treatment of the fermionic case is also 
illustrated by the cancellation of logarithmic singularities of the iterated ladder and stretched 



box diagrams of the Yukawa model 
of the box diagram 32|. 
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32l | . This cancellation yields the full covariant form 



VIII. CONCLUSION 

In this paper we propose a conserved electromagnetic current operator that acts on the 
valence component of the three-dimensional LF wave function of a two-fermion system. 
In order to obtain the LF current, we have exploited the quasi-potential approach to the 
BS equation [8] and then we have projected the relevant quantities onto the light-front 
hyperplane. This approach has been already applied to the BS equation for both boson and 
fermions and to the construction of the conserved current operator for two-boson 

systems [16]. 

The starting point of the reduction scheme to the LF is the QP approach to the BS 
equation with a proper auxiliary four- dimensional Green's function without instantaneous 
terms. Such terms, that represent the peculiar feature of a fermionic system, are recovered 
both in the kernel of the LF equation for the valence wave function and in the LF two-body 
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current operator. The integration of the minus component of the particle four-momenta in 
the relevant operators, namely the Green's function, the T-matrix and etc., allows one to ac- 
complish the desired LF projection of both eigenequation and em four-current. Furthermore, 
we have identified a non-unitary operator, the reverse LF-time projection operator, that acts 
on the LF valence wave function and allows one to reconstruct the four-dimensional BS 
amplitude for bound or scattering states. Such an operator is non unitary since the valence 
component does not carry the full normalization of the wave function (for a discussion of 
the probability of the higher Fock-components in the Wick-Cutkosky model see Ref. jl?!). 
The explicit expression of the reverse LF-time operator in terms of the quasi-potential is 
essential for obtaining a LF three-dimensional current that fulfills the WTI, at any order 
of the quasi-potential expansion. In particular, the reverse LF time operator applied to 
a four- dimensional current generates a three-dimensional LF current. Such an application 
trivially leads to the equality between the matrix elements of the four-current, evaluated by 
using the full BS amplitudes, and the matrix elements, of the LF current evaluated by using 
the valence wave functions corresponding to the previous BS amplitudes. Then the current 
conservation follows. This result is essential for identifying the ingredients to be used for 
demonstrating the WTI for the truncated LF current. In particular, we have defined the 
left and right LF charge operators (that do not contain interaction), and we have found the 
formal expression for the LF WTI, where the LF Green's functions, for the initial and final 
states, appear. The last step in our analysis, is of particular relevance. Indeed, a naive 
truncation of the QP expansion that defines the LF current doesn't lead to a conserved 
current, while retaining all the contributions up to a given order in the effective interaction 
allows one to construct a LF truncated current that satisfies the WTI. In such a truncated 
WTI, the truncated initial and final Green's functions appear, as well as the same LF charge 
operators obtained in the full case. In particular, the truncated effective current operator 



can be put in correspondence to a sum over intermediate states [Ul] up to some maximal 
number of particles exchanged at a given LF time that flows from the initial to the final 
valence states due to the photo-absorption process. The truncation of the expansion of the 
quasi-potential implies that the observables derived from the three-dimensional conserved 
current operator are frame dependent, but with gauge invariance correctly implemented. 
The covariance under Lorentz transformations of the full LF current operator, j'^{Kf,Ki), 
is not greatly relevant, given the equality in Eq. (1571) between the four-dimensional matrix 
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elements of the covariant current and the corresponding three-dimensional ones. However, 
since the truncated theory is necessary for obtaining workable approximations, the covari- 
ance properties of both the full and the truncated LF current will be investigated elsewhere 



251]. Here, we can anticipate that the covariance under the seven kinematical LF transfor- 



mations is satisfied by the LF current operator (and also by the valence wave functions) in 
both cases, with a suitable introduction of new factors in the vertical bar operation. For the 
truncated theory, the violation of covariance under dynamical LF transformations produces 
effects on the matrix elements of the current, that can be reduced at any desired accuracy, by 
increasing the order n of the truncated quasi-potential, namely approaching the full theory. 
This necessarily leads to consider intermediate Fock states with larger and larger number of 
particles in the evaluation of the effective interaction. Moreover, a quantitative study of the 
effect of covariance violations was already performed in the computation of the masses of 
bound states for the bosonic model j^, [2^ • It was concluded that the expansion in the Fock 
space is rapidly converging for a given covariant model (see also j^, 0, 30]). In particular 



the splitting between magnetic states of spin-1 composite bosons decreases when the kernel 
of the LF bound state equation is evaluated taking into account higher order terms, like 



stretched boxes, in a bosonic model 
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Our procedure for constructing a truncated LF current has been illustrated in an actual 
case: the Yukawa model with chargeless boson exchange in ladder approximation. We have 
evaluated the LF current operator at the lowest nontrivial order and have explicitly checked 
the Ward-Takahashi identity for such model. The role of instantaneous terms has been 
clarified and their relevance in producing two-body contributions has been emphasized. 

Let us finally comment on possible problems about singularities and regularizations that 
can occur in the truncated LF current. Restricting to a model without self-energies and 
vertex corrections, but considering in the kernel the ladder and two-body irreducible cross- 
ladder terms, the Bethe-Salpeter equation is finite allowing a solution. However, the proiec- 



tion of this equation to the light-front is plagued by infinities (see, e.g., [lO|, 1291, [3^, ISll. l32|). 
Recently it was shown that the finite covariant box amplitude in the Yukawa model is fully 
recovered when all terms in the light-front projection beyond the iterated box (i.e., stretched 
box and instantaneous terms) are obtained and that all the singularities are canceled (321] . 
Thus in general one could argue that regularization problems occur in the integration loops 
for the truncated current at order > 1, since the presence in Aq of the global propagation, 
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described by Go{K), destroys the balance between singularities that allows for a finite result 
when only Gq{K) is considered (see e.g. 32|]). For instance, an easy form of regularization 
can be introduced in Go{K) through a cutoff function, e.g. 6'(yU^ — Mq). Since the starting 
four-dimensional covariant BS model is finite, one expects that the effect of the scale n 
vanishes if the QP expansion is not truncated. 

However, the issue of renormalization of the nonperturbative bound state problem in 
the truncated QP expansion is subtle |2J|, and therefore the dependence upon the scale fi 



in the three-dimensional truncated theory should be carefully analyzed. This non trivial 
elaboration on the renormalization issue has to be postponed to a future work. 

In summary, we have proposed a systematic expansion of a conserved electromagnetic 
current operator within LF dynamics for two-fermion interacting systems, using the quasi- 
potential approach to the Bethe-Salpeter equation. As to the future perspectives, we plan to 
apply such an interacting current for the investigation of inclusive and exclusive electromag- 
netic processes, like hadron form factors and deeply virtual photon scattering, after properly 
generalizing the present approach to fermion-antifermion systems 25|]. It should be pointed 
out that the identification of the matrix elements of any operator in the four- dimensional 
space with the matrix elements of the corresponding three-dimensional operator acting on 
the valence wave function is a general procedure. 



Acknowledgment s 

This work was partially supported by the Brazilian agencies CNPq and FAPESP and by 
Ministero della Ricerca Scientifica e Tecnologica. J. A. O. M. and T. F. acknowledge the 
hospitality of the Dipartimento di Fisica, Universita di Roma "Tor Vergata" and of Istituto 
Nazionale di Fisica Nucleare, Sezione Tor Vergata and Sezione Roma I. 



APPENDIX A: USEFUL IDENTITIES 



We introduce the following identities that will be useful in exploring the relation between 
LF and four-dimensional covariant quantities. It is straightforward to get that: 

(l^ - m)(^„ + m) = (p^ - m^) — (|^o„ + m) , (Al) 

Zp 



29 



{ion + m){i-m) = {p^ - m'^){j>on + 'm)—— . (A2) 
These identities imply the following relations 

{k[-\G,\K)GoiK)\k^) = ^^(^,^^ + m,)i^2on + m,)5ik[- -ki) , (A3) 

{k[-\G,{K)G^\K)\k^) = {hon + mi)(|^2o„ + m2)^^6{k[- - k^) . (A4) 

Note the presence of the product 7j*'7^ that is essential for the definition of the LF charge 
operators, Eqs. fl5^ and dHSD- Exploiting (7"*")^ = 0, one obtains another useful identity 

^+ 

(l^ + m) — (i^ + m) = (^„ + m), (A5) 
Zp 

with = (p^ +m)/p+. 

Using i) Eqs. (lASP and (lASp . and ii) the explicit form of Go(-^), one gets 

g,{K) = \G,{K)G^\K)Go{K)\ , (A6) 

for the free three-dimensional propagator. 

Now we can relate the interacting LF Green's function directly with the four- dimensional 
Green's function by evaluating 

\G,{K)G^\K)G{K)G^\K)G^{K)\ = \G^{K)G^\K)Go{K)\ + \G^{K)T{K)Go{K)\ 

= go{K)+go{K)Uo{K)T{K)Uo{K)go{K) 
= g,{K) + g,{K)t{K)g,{K) = g{K) , (A7) 

where Eq. for the definition of t{K) was used. 



APPENDIX B: THE INTERACTING REVERSE LF PROJECTION OPERATOR 

In this Appendix, we will prove some useful identities involving the interacting reverse 
LF projection operator. 

The following relation between n(i^) (see the definition in Eq. ( H2i) ) and T{K) can be 
obtained from Eqs. (EI]), ([27]), ([39]), ([29]) and (El) Indeed one has 

U{K) = Uo{K) + Ao{K)W{K)IloiK) = 

= Uo{K) + Go{K)W{K)Uo{K) - Uo{K) go{K)Uo{K)W {K)Uo{K) = 
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Uo{K) + GoiK)W{K)Uo{K) - Uo{K) go{K)w{K) 
Ii^{K) g,{K)g-\K)+G^{K)W{K)Il,{K) g{K)g-\K) = 

[Il,{K) + G^{K)W{K)Il,{K) + G,{K)W{K)Il,{K)g^{K)t{K^ g,{K) g-\K) = 
1 + G^{K)W{K) + Gq{K)W{K)G^{K)T{K)\ Yi^iK) go{K) g-\K) = 
[1 + Go{K)T{K)] no(ir) go{K) g-\K) (Bl) 



Furthermore, by using Eqs. ([2]) and f[T9l) . one gets 

Ii{K) = G{K) G,\K) G,{K)\ g-\K) (B2) 
that can be recast in the following form 

G-\K) Il{K) = G,\K) Go{K)\ g~\K) (B3) 
The analogous expression for Il{K) reads 

n{K) G-\K) = g-\K) \G,{K)G^\K) (B4) 

APPENDIX C: ZERO-ORDER LF CURRENT OPERATOR AND WTI 

In this Appendix we will evaluate the matrix elements of the free LF current operator, 
Eq. (!69l) . between free particle states and we will prove explicitly the WTI. The matrix 
elements of the free 4-dimensional current operator are 

{ki\J^mPi) = -27r [ei 7^ 5\h - Pi - Q) ((/^/ - h) ' "^2)] 
^\\^2,k^^Ki-k^,Vl^i^^-Vl\ , (CI) 

where = — . The factor (— 27r) is introduced in the current operator to make it 
compatible with the free Green's function, see Eq. ([3]). 

Since g^^ is the identity in the two-particle space, modulo some factors, to simplify the 
presentation let us consider in what follows only the relevant part of j'^^^^'^'> containing the 
^"-integration. Moreover, in order to make more fast the discussion related to the position 
of the poles, we will take profit of the 6'-functions present in the dropped g^^ factors, that 
remind us we are dealing with particles in the external legs. Therefore, one has 

e {k'+k[^\ \GoJ^Go\ \kthi.) = 
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1 r e,5 {k[+ ~kt- g+) 5^ (k[^ - h± - Q±) 



(2vr) 



' ' k[^kt {k^ + Q-- k[-^ + z^) (fcr - fci", + 



(^U + ^i)7f (^io« + - + ma) , ^ , , o (C2) 

where fc; = ki + Q, and = ^(A;+) 9{kl')e{K^ - A;+) - k^'). It is important noting 

that, without profiting of the presence of 6 on the left of Eq. ( ]C2p . a lengthy discussion 
of the poles leads to the presence of 6 in the result. We have used the identities (lAip and 
f lASP to simplify the following combination that appears in the numerator of Eq. (]C2p 

ihon + fn2)i^2 - m2)i^2on + ^2) = (kj - ml){^2on + ^2) , (C3) 

with k2 = Ki — ki = Kf — k'^. 

Integrating over k^ and assuming that > and > 0, without loss of generality, 
one gets that 

e {k'+k[^ I \GoJo^Go\ I k+h^) = iQ e,5 {k'+ - k+ - Q+) 5^ (k[^ - h± - Q. 

X 



k[^{Kt-kt)ktAf-Af^^- 
+ 1^2 (C4) 

where we have introduced the four- vector quantities A^^-* = Kf — k[^^ — {Kf — k[)on + 
and Af^ = Ki — kion — {Ki — A;i)on + for convenience. Note that the minus component 
of A/ and Aj are the only non vanishing ones. The on- minus-shell values of the individual 
momenta are 

'^lon — ^/+ \^'~') 

(K h'\- - Aa.? + ^1 .p7^ 
(^/ - ki)on = _ (^7) 



(i^.-fci)o-n= ^^"\ 7"^^ - (C8) 



{Kj ^ - + ml 

Kt - kt 

Taking into account the definition of go{K), Eq. (fT6ll . and the matrix elements of ei^p, 
Eq. fl65|) . we have 



^7712 

+1^2, (C9) 
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Finally, by multiplying by the proper go{K) ^, the matrix element of the free current oper- 
ator can be written as 



(K^K^ \f^'^'\ Kkii.) = -ieK{k[^MK{hon) 
Kt-kt 



2m2 



-{k[^,k[^\ei,LF\ktkii.) + 1^2 



(CIO) 



To derive the WTI for the free LF current operator, we evaluate the four-divergence of 
the current by contracting it with Q^. Since {Kf — k[)~^ = {Ki — ki)~^, which follows from 
kinematical momentum conservation, the momentum transfer can be written in terms of 



A's as, Q = Af^ — A-°^ + k[^ — kion, and thus 



2_ (^(0)- _ ^(.)- 



(Cll) 



Then one has 



A+(fcU)^A+(fcion) = A+(fcU)^A+(fci,„,)(Af- - Af-) . 
/^From the above results, the four- divergence of the free current becomes: 



(C12) 



{k'+k' 



i± 



Q-j 



;{0) 



X 



kth±) = -te A+{klJ^A+{k,on)A+i{K, - k,)on) 



Kt -k t,^ (0)- 
2m2 

{k'+k[^ [goiKfT'Q^.^-QfMm^' ktki±) 



{A';>- - Af'>-){k[+, k[^\ei,LF\ktki,_) + 1^2. 



(C13) 



since the free resolvent that appears in Eq. (]C13|) comes from the following relations 



Af)'A+(A;i„„)A+(A;2,„)=2 



2mi 2m2 



-1 



Af -A+(A;U)A+(fc2on) = ^^ 



kt K - kt 
2mi 2m2 



Equation flC13P gives, as expected, the matrix elements of fITOl) . i.e. the WTI for n = 



(C14) 



APPENDIX D: WTI FOR THE FIRST-ORDER LF CURRENT OPERATOR 



In this Appendix we show the WTI for LF current, obtained by truncating at the first 
order the effective interaction, see Sect. VI. Let us rewrite the LF first-order current operator 



:di) 
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Using the current conservation for Jq and JT/, see Eqs. ( ITT!) and (fT2!) . the four- divergence 
is given by 

Q'^jf ^ = Qo'^QIf - Qlf9o~^ + no [{eV - Ve) + VAoGQ-^e - eG.'AoV] Uo . (D2) 

Note that the term [—V Aq e Gq ^ Ilo + no Gq ^ ^] is not present in the above equation, 
since it is vanishing due to the absence of a Dirac structure in the operator fT and because 
of Eqs. (jMD, (jA4]) and (15^. 

Furthermore, by using the definitions of i) Aq, Eq. ([21]), ii) Gq, Eq. (1271) . and iii) 2^;^, 
Eq. ([62]), and Qf^, Eq. ([63]), one has 



M.-c(l) 



V Go Go^e-eGn^ Gn 



'-^o 



^7o"'QiF - Qlf9o-' - no [\/ Ho go Hq Gq - cGq ^ Hq (7o TTq Hq 
go-'Q^F - Qlf9o'' - [w^'^ \Go Go^ e Ho - Ho e Gq^ GqI w^'^] = 



w 



(D3) 



APPENDIX E: WTI FOR ARBITRARY n > 1 LF CURRENT OPERATOR 



The proof of WTI is based on the induction hypothesis. Once the current operator at 
order n satisfies the Ward-Takahashi identity, we must demonstrate its vahdity for n + 1. 
First note that 



J 



WnAoJj^ + Ji'AoWn + Wn+lAoJo^ + Jo^ AoW^+i 



n-1 



(El) 



1=1 i=l 

where we have separated out the free term from the interacting one in the 4-dimensional 
current using Eq. Since Q^Ji = eV — Ve, by induction, we have: 



g-'QiF-QLF9-' + ^o 



WnAoieV - Ve) + 



{eV - Ve)AoW^ + W^+,Ao{Go'e - eGo') + (Gq - eGo-^)Aoiy„+i + 



n-l 



W,Ao{eV - Ve)AoWn-i + ^^^o(Go - cGq ')Aoiy, 



i=l 



i=l 



Ho . (E2) 
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In the above equation, the term Ho [— iy„+i AqcGq ^ + Gq ^eAoWn+i] Hq vanishes due to the 
absence of a Dirac structure in the operator e and because of Eqs. (IA3p . flA4p and fl52]) . By 
using the same relations and the definition of Go, Eq- (l27j) . we can simphfy the following 
term that appears in last sum, i.e. 



Ao(Go - eGo ')Ao = (eAo - Aqc) 



(E3) 



Then, since by definition l/AoWj = Wj+i, we have 



V J/t — Un ^LF ^LF iJn + ■'-'-0 



n-l 

ly^+iGoGo + eGo 'Gol^n+i + ^i^oeWn-^+l - W,+leAoWn-^ + 

n 

-i+l 



n-l 



i=l 



i=l 



i=l 



(E4) 



and rearranging the terms we get 



(-)^i-c{n+l) _ -1 qL _ qR -1 I TT 



n-l 



eGo^Gol^n+i - Wn+iG^G^^e+ 



i=l i=0 

Since the sums cancel each other, we have 

^l ■c(n+l) 



i=l 



(E5) 



, 9n' QIf - QIf 9n' + no [eGo ^GoW^n+i - l^„+iGoGo-^e 

gn' Qip - Qfp 9n^ + QlfWu+I ' Wn+lQ^p 



Tin 



(E6) 



where the LP charge operators as given by Eqs. (l63l) and (162|) have been introduced through 

noeGo-^Gol^n+l Ho = Q^pWn+l 

UoWn+iGoG^'eUo = Wn+iQ^F , (E7) 
Finally, since g:;^^^ = g-^ - one gets 

,M,c(n+l) ^ g-l^ qL^ _ qR^ _ ^E8) 



Thus, by induction, we conclude that the LP electromagnetic current operator j^*-"^ is con- 
served at any given order n of the quasi-potential expansion, once the matrix elements are 
taken between eigenstates of: fif^^l^n) = 0. 
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APPENDIX F: INTERACTION-DEPENDENT PART OF THE LF FIRST- 
ORDER CURRENT OPERATOR IN LADDER APPROXIMATION 

The contribution of the interaction to the LF current operator in lowest order in ladder 
approximation comes from two-body irreducible amplitudes given by the second term in the 
r.h.s of Eq. flHT]) . Using Eq. fl3T]) for Aq one has: 

where w^^^ = IIo ^ Hq is the three-dimensional effective interaction (see Eq. fl83|) ). 

In the following we analyze the matrix elements of the first-order current operator relevant 
for the em processes in the spacelike region. 



1. Evaluation of the term \GoJ^iQ)GoVGo\ 

Let us start with the third term in the r.h.s. of Eq. (IFip . gQ^\GoJ'Q{Q)GoVGo\ g^^, 
which is diagrammatically illustrated in Fig. 2. It is important to observe that in this 
term contains the relevant pair production contribution, for > (cf. diagram (b) in Fig. 
2). For the sake of simplicity, we can first drop the multiplicative factors Qq^ on the left 
and on the right, and we will consider the current only for particle 1. However, we can take 
advantage of the ^-functions contained in both g^^{Kf) and g^^^Ki), that greatly help in 
the discussion of the analytical integrations over and k'f . Therefore, we have to evaluate 
the following matrix element between free particle states 

&{k'+k[^\ |GoJo''(l)G'ol^Go| \kth±)=^el(^^^ 6 



X 7f 



(Pl - Plan + iK - Pt) -Pi - (^1 - Pl)on + ^I+Z^) 

l)on 

^^^^^ - {K - Kon + - K - {K, - h)on + ^J^) 



where 6 = e{k+')e{K+ - A;+) e{K^ - kf ) and p'^ = k'^ - Qf". Eq. ^ is represented 
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by the Feynman diagram shown in Fig. [TJ The six poles in Eq. (lF2p are 

Kb = - {Ki -ki)-^ + i ^ 



Kt-kt 

e 



kj -Pi 



Pi = Kc -Q' = Plan - (F3) 

Pi 

with the on-minus-shell definition of the respective momenta given by 

kjj_ + ml 



{K, - k, 
[ki-pi] 



kt -Pi 

p^j_ + m( 



Plan + 
Pi 

, ,„ _ k[\ + ml 



= (F4) 



The integrations over fcf and k'^ in Eq. flF2l) are performed analytically using Cauchy's 
theorem with the conditions Kt > and > 0. As a result, we get the six contributions 
that appear in Fig. [2l Note that diagrams (d) and (e) are two-body reducible terms which 
are canceled by us now discuss in detail the diagrams of Fig. [21 

devoting a specific subsection to the pair diagram (b). 



a. Diagrams (a), (c), (d), (e) and (f) for = k't' — > 

Let us start considering the kinematical region p^*^ = k't — Q^ > and kt > kt — Q^ ■ We 
first perform the analytical integration of Eq. (IF2p on k^ . The poles /c^ and k^(j belong to 
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the lower semi-plane, while k^^ lays in the upper semi-plane. The result obtained integrating 
in the upper semi-plane contains contributions corresponding to diagrams (a), (c) and (d), 
viz 



Q {k[^k[j^\ \GoJ^{l)GoVGo\ |A;+fci±)(a)+(c)+(d) = ei 
1 + mi 



e {pt) e {kt - pt) 



X 



X 



dk 



271 



1 1J+ 



-I- i-^ 

Ion ^ fc^ 



Kj - k[- - {Kf - k[)~^ + tj^j^ 



{kt - pt)Pi 



Kj-k[--{K,-k,)-^ + z 



X 



X 



^lon + mi 



[k[- -Q-- Pi, 



e \ 2 

Pi 



(^1 -Pi)on + «r+ 



kt{Kt-kt) 



(^io„ + mi) [{fCi - + ma] 



(F5) 



{Ki ki),^ + «^+^^+ ki,^ + 

Note that i) the instantaneous term, proportional to 7^*", leading to a first contribution to the 
current (c) shown in Fig. [21 is explicitly separated out in Eq. f IFSp . ii) the global propagation 
of the initial state, i.e. 1/ [K~ — {Ki — ki)^^ — fcfo^] will be canceled by go{Ki)~^. 

In order to separate the processes corresponding to diagrams (a) and (d) of Fig. [21 we 
make use of the identity 

1 1 

{kj - k'f - {K, - ki)-^ + - {ki - pi)-^ + (k[~ -Q-- p^on + 



+ 



(^Kj -k[ - {Ki - ki)^ + i-^^fzq; - (ki - + i 



K^ -Pion + i^ - {Ki-ki)-^ + i-^ 

Pi '^i 



(^1 - Pl)on + i 



[k'- -Q- -p-^^ + t^Y 

(F6) 



K~Pi 



and integrate Eq. (IFSP analytically on k'^ . The first term in the square brackets in Eq. (]F6P 
generates the contribution illustrated by diagram (a) in Fig. ([2]), once we take the residue at 
the pole k'lj^, in the lower semi-plane, and multiply on the left by gQ{Kf)^^ and on the right 
by go{Ki)~^ (cf Eq. ( TT71) ). The second term in the square brackets leads to the contribution 
of diagram (d), a two-body reducible term, which is canceled out by one of the contributions 



in j'^^'^^^gQW^^\ see Eq. (IFip . The diagram (a) is given by 



{k'+k[^\j<'^^\ktki^)^,^ = zei{2mi){tgf 



9 {kt - pt) (pt) 
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(FT) 



The four-vectors A^}""^ and A-'^"'' are the following combinations of the four-momenta 



(la) 



A 



(la) 



- Kon - i^i - kl)on - ih - Pi)on + iS 
Ki - Plan - {Ki - ki)on " (^1 " Pl)on + is . 



(F8) 



The minus components yield the three-body (2 fermions and the exchanged boson) global 
propagation in the final and initial states, respectively. 

The instantaneous term (c), in the kinematical region under consideration in this subsec- 
tion, can be obtained from the same pole A;'^^. Then one has 



Pi 



1 on 



(F9) 
We first 



Let us then consider the kinematical region where = k'^ — > kf. 
perform the analytical integration of Eq. (lF2p in k^ by calculating the residue at the pole 
in the lower semi-plane, where only this pole appears. Then in this region we obtain two 
contributions from Eq. (]F2|) : a reducible one (cf diagram (e) in Fig. ([2])) and an irreducible 
one given by diagram (f), namely an instantaneous term. Diagram (e) is canceled by the 
remaining part of j^^^^^gQW^^\ The result corresponding to diagram (f) is 

e {k'+k[^\ \GoJo^iQ)GoVGo\ \kth±)if) = ei^ 6 9 (pf) 6 {pf - kt) 



X 



X 



X 



^lon + "^1 



(/T/ - lt[)on + m2 



Kj-k[--iKj~k[)~, + t 



1 



7r 



ipt - kt)vt \k[- -Q-- k-,^ - {p, - fcOon + + 

it ^l^2a {hon + mi){{fii-^i)on + m2) 



2 kt{Kt - kt 



^lon + ^1 



(FIO) 



The analytical integration of Eq. flFlOp is easily performed by calculating the residue at 
the pole k'^Q in the upper semi-plane. Therefore, the contribution from the instantaneous 
term for k't' — > kt, represented by diagram (f), is: 



(k'+k' h-W'^l/t+fc. , Vn - -e.Ua? e ^ - 

(fci fci^lj \k^ki^)i^f)- ^ei{ig) U {p+-kt)A^^^' 

X AAKonH7t^l^+{klon)A+{iKf - k[U)T2^A^m - h)on) 



(Fll) 
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where in the denominator appears the minus component of the following four-momentum 

Af ^) = Ki- hon - {Kf - k[)on - {Pl - h)on (F12) 

which yields the intermediate propagation of the three-body system composed by two 
fermions and a boson. 



b. Diagrams (b) and (c) for > — Q = Pi 

In the kinematical region where > k'{^ — = pf the non vanishing result can be 
obtained from the pole k^^ in the upper semi-plane, as in Eq. (IFSp . Then one has to 
consider the pole k'{^, when the integration over k'{~ is performed. Note that in Eq. (IFSP 
0{pi) has now to be substituted by 6{—pi). The resulting matrix element of the current 
operator (b) is 

Pl K"^! Pl I f 7 

X A+(A;U)7i'A+(pi„„) r-A+(A;i„„)A+((ir; - A;i)on)r2„A+((K, - fc^on) , (F13) 

where the combination A^^*"* of four momenta is 

Af) = Q-A;U+pi,„ + z£ . (F14) 

Finally one has to compute the two-body current due to the instantaneous term (c), 
which in this kinematical region is given by: 

X A+iKonhi it T'^A+{han)A+{{Kf " k[U)T2aA4K^ - k,U) . (F15) 



2. Evaluation of the term \GoVGoJ^{Q)Go\ 

Let us now discuss the first contribution in Eq. (IFip . go^\GoVGoJ'Q{Q)Go\gQ^, depicted 
by the diagrams (g) to (i) in Fig. [3l It is worth noting that the pair contribution is not 
present due to the conservation of the plus momentum component. For the sake of simplicity, 
the two reducible diagrams, which are present in this first term and are canceled by the term 
w^^^goj''^^^'^ of Eq. (IFip in an analogous way as we have already discussed for diagrams (e) and 
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(d) of Fig. [21 are not shown in Fig. ([3]). As we did for the third term in Eq. (IFip . we evaluate 
first the matrix elements of the relevant part of the contribution under consideration, i.e. 

1 Kon + ^1 (/T/ - ft[)on + 1712 

riFaa fi + mi ^ 

(A;'+ - {k[- - p[~ - {k[ - p[)-^ + z^^) {p[- - p[-,, + ^4 ' 

J- f^lon \)on 

^^'■(^^ - (^r - Kon + ) [kt - k-, - {K, - + 



where the intermediate momentum is changed now to p'^ = k'^ + Q'^. In order to analytically 
evaluate the integral, in addition to the poles /c^"^, k^^, k'{^, k'^^, given in Eq. (1F3[) . we have 
to consider the following ones 

^ID + Q = P'lon ~ "^"Tf 

Pi 

Kd = P'r + iK - p'i)on - ^77+ 7+ (F17) 

- Pi 

with the on-minus-shell definitions given by 

_ + mj 

Plan /+ 
Pi 

— 

As already said in Sect. VII, we have always p'^ > and therefore no pair contribution can be 
generated from this term. In what follows we consider two regions: i) k'{^ > p'{^ = kf + 
and ii) p'^ = kf + > k'^ . Moreover the definition of Pi suggests to start with the 
integration over k'^ . In this way we have to discuss only three poles. 

a. Diagrams (g) and (h) for k'-^ > p'^ > 
From the residue at the pole k'{^ one gets: 

e {k'+k[^\ \GoVGoX{l)Go\ \kth±)(9)+ih) = & e{k'^ - p'+) 
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f^lon + '^l 



[if<f-lt[)on + m2\ 



X 



2a 



- i^f - K)on - P'l - iK - K)on + ^^fZ^ + ^ K+-k'+ 



X 



X 



fion + mi 



+ 



Pi J 
Ion + mi 



7i 



^^^^^ - (^r - Kon + ) (^r - K - {K. - A:i).-„ + 



(F19) 



Adopting the same strategy as in Eq. (\F6\i . one obtains i) the contribution illustrated by 
diagram (g) in Fig. ([3]) and a reducible term (not shown in Fig. ([3]) ), ii) the instantaneous 
term (h). In order to separate the irreducible contribution from the reducible one, we exploit 
the following identity 

1 1 



KJ - (Kf - k[). 



p'l 



{k[ — p'i)on + ^I7T~7T 



Pi ~ Plon ' 



i-^f ^l)on Pi P'l)on + ^fc^+^p'+ + 



Pi Plon + ^ jj+ 



1 



X 



KJ - {Kf - k[ 



'on -Pi on 



p'l) on + ^A.'+^„'+ + ^ 



KJ-kj^ Pi 



(F20) 



In particular, the term generated by the first contribution in the square brackets can inte- 
grated over kj taking the residue at the pole kj^^. Then one gets 

e{k[^~p'i^) 



{k[+k[^\f^'^'^\ktki^)^,.=zei{2mi){zgf 6 



xA+(A;U)r?A+(Kon)7rA+(A;i„„)A+((K^ - k',),^)r,^A^iiK, - kiU) , (F21) 



where the combinations A^^^^ and Aj^^^'' of four momenta are 



A^;^) = Kf- pI^ - {Kf - k[U - (K - p'i)on + te 



,{19) 



A^^^^ =K,- ki,r. - {Kf - k[),r. - {k'l - p'i)on + ^£ = A 



(1/) 



(F22) 



since k'-^ — p\ = pi — ki. The second term in the square brackets yields a reducible contri- 
bution, which is canceled out by a term from w^^^goj^^^^^. 
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The result for the instantaneous diagram (h) can be obtained by taking the residue at 
the pole k^j^ in Eq. fIFlQp . i.e. 

9 {k[^ - 



{k[^k[^\j<'^^\kth±)iH) = ^e,{zgf e ^^^^ 



ir2,A+((K, - k^) 



(F23) 



h. Diagram (i) for p'^ > k'^ > 

In this kinematical region only an instantaneous contribution is an irreducible term, the 
one illustrated by diagram (i) in Fig. ([3]). By taking first the residue at the pole k'{~^ and 
then the residue at the pole k^^, one has 

^ [Pi - ^1 )Pi A) 

xA+(A;U)r?7i+7rA+(^ion)A+((i^/ - A;;),„)r2„A+((K, - k,),^) , (F24) 

where the three-body denominator is the minus component of 

A^^^) = Kf- kl^ - {K, - k,U - (P'l - K)on = Af^^ , (F25) 



APPENDIX G: CURRENT CONSERVATION IN THE LADDER YUKAWA 
MODEL 

The operatorial form of the current conservation at the first order of the QP expansion, 
as given by the expression in Eq. (!74l) . will be explicitly checked in what follows. The 
symbol = 9{ki)9{k'^)9{Kf — ki)9{K^ — k'^) simplifies the notation. Furthermore, it 
is worth noting that the A's, given in Eqs. ([FS]), flFT2D . flFTlD . flF22D and (IF25D . allow to 
express the momentum transfer, Q, in a convenient way, depending upon the term we are 
considering, and they have only minus components non vanishing. 

Let us separate two kinematical regions : i) k'^ > and ii) > k'^. 



1. p+ = k'+ -Q+ >0 

In this region one has contribution from diagrams (a), (c^), (f), (g), (h) in Figs. ([2]) and 
[l3]) (since k'^ — p'l = Pi — kf > 0). For diagrams (a) and (g), the momentum transfer can 
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be expressed in two different ways, viz 



— lAr LA^ -t- pion , 



Q = ^f'^-^r-hon+pln- 

Therefore, since A+{k[^J{lt[^^ - ^ion)A+{pion) = 0, one has 

A^{klJQA4p,on) = A+(A:U)|^A+(pion)(A(^'^)- 
Analogously one finds that. 



(19) 



-A, 



(la)- 



)• 



A, 



(19)^ 



(Gl) 
(G2) 

(G3) 

(G4) 



Eqs. (1G3P and (lG4p allow to calculate the four- divergence of the current. 

The divergence for the term (a) plus (g) (see Eqs. flF7p and flF2ip ) is given by: 



{kt - Pt) 
9 - / 1 



XA+{k2^jT2aA+{k2on) 



-(ip ) A+(A;U)^7i+A+(Pio„)r?A+(A;i,„) 
A}^ J Pi 



A+(A;U)r?A+(p;„J^7+A+(A;i,„) 
Pi 



A 



(G5) 



where we have used /c2on = {Ki — ^i)on and /c2o„ = {Kj — k[)on to simplify the notation. 

Let us consider the divergence of currents (c^) and (i^), corresponding to the contribution 
of instantaneous terms in the kinematical region k'-^^ > Q'^ > 0. The following relation, 
obtained by exploiting i) the property (7"*")^ = and ii) the anticommutation rule {^, 7"''} = 
2A;+, is useful for the next formal steps, i.e. 



^Ak'ionm^ = K{k[J - ii] 7+ = A^{k[J [^l, - i>,J, 7+ 
= A+(A;'i„„) [mi -i^ion] 7^ = A+(A;io„,) [-2p| + 2mi7+A+(pio„)] 
= 2p+ A+ 



-1 + — 7+A+(pi„„) 
Pi 



(G6) 



Then, from Eq. (lF9p . one has 

X A+(A;^Q„ 



fi[k't -Q^)e{kX-vX) 



la- 



-1 + — 7+A+(pi, 
Pi 



(kt-pt)^} 

T'^A+{k,on)A+ik'^jT2aA+{k2on) ■ 



(G7) 
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For diagram (i), Eq. (lF24p . in this region, a relation analogous to Eq. (1G6|) holds, i.e. 



7+ Q A+(fclon) = 7 [A - h] A+(A;ion) = 7 \Aon " hon] A+(A;ion) 



2pi 



1 - ^A+(p'ion)7^ 



(G8) 



Therefore, from Eq. (]F24p . one obtains 



Q-{k[+k[^\j<'^\kth^)'=iee,{ig) 



{kt-pt)Af^^- 



1 - A+ Pio„ — 7i 
Pi 



xA+(A:U)r? 
Summing Eqs. (IG7p and (IG9p . we find that 



Q.{k[+k[^\J<'^\kth±)U^^.)=^ee^{^g) 



A+{kion)A+{k2on)'^2a-^+{k2on) 



,e {k[^ - Q+) 9 {kt - pt) 



(G9) 



xA+{k[ 



Ion) 



Pi Pi 



{kt-pt)Af^^- 

A+{kion)A+{k'2on)'^2aA+{k2o4G10) 



The instantaneous terms (f) plus (h), see Eqs. (IFlip and (]F23P can be calculated following 
the same steps as the ones that lead to Eq. fiGTOl) . Noting that Af^^ = A[^^\ one has: 

,e {k[^ - Q+) e {pi - kt) 



Q ■ {k'^k[^\j<'^\ktki^)^f)+^,) = I e e,{ig) 
xAj^ik^g^ 



^7i+A+(pi.„)r? - r?A4p;,j^7 

LPi Pi 



ipt - fc^)Af 

A+ ( A;io„) A+ ( A;^„„)r2„A+ ( A;2on()Gl 1) 



Adding Eqs. (IGSp . (IGIOP and (1G11|) one gets in the kinematical region k'^^ > > the 
result 

Q ■ {K^k[^\f^'^\kth±) e {k[+ - Q^) = 



e{k'+-Q+) e 



ei A+ ( A;;, J ^7i+ A+ (pio„) {ptpiM'^ ^ \ kth±) 
Pi 

ITLl 



- {k[^k[^ 1^(1) (Kj) biVix) ei A+ {pU -liltAAhon) 

Pi 



(G12) 



where the matrix elements of the effective interaction in first order of QP truncation, for 
total momentum Ki, are given by: 

{ptPi±\w^'\mkth^) = z{zgf e{pt) e{kt)e{Kt - pt) e{Kt - kt) 

^{pt-kt) 1 , 9{kt-pt) 1 



ipt-K) A^'^- (K-pt) A^)- 



A+(pio„)r°A+(fcio„)A+(4„„)r2aA+(A;2o„) 

(G13) 
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For total momentum Kf^ the matrix elements of w^^\Kf) can be easily obtained from the 
previous expression properly changing the individual momenta, i.e. 

{k[^k[A^^'\Kf)\p\^A^) = liigf 9ip[^) 9ik[^) e{Kj-p\+) 9iKl - k[^) 

0{k'+-p[+) 1 e{p'+-k'+) 1 



X 



^+{k[on)'^l^+{Plon)^+{k2on)'^2aK{k2on) 

(G14) 



Let us remind that k'^ —p'l = Pi —k^ and 9{p^) is redundant. For Fq, = 1, i.e. for a scalar 



boson exchange, one obtains the expression derived in Ref. 



io|. 



2. Q+ - k'+ > 

In this kinematical region diagrams (b) and (c^^) contribute. Furthermore we consider the 
contribution of diagram (i) in this region (cf the discussion after Eq. (1F24P ). The divergence 
of the pair current (see Eq. (IFlSp ) is given by: 

,9 {p[^ - k[^) 9 [Q^ - k[^) 



Q ■ {k[+k[^\j<'^\ktki^)^,) = -z e eiiigf 



xA+{k[^j'^^t^+{Plon)^A+{kion)A+{k2on)^2aA+{k2on) , 
Pj 



(G15) 



where Q = A^^^ + k[^^ — pion has been used (see Eq. (IF14p and note that only A^^''-'" does 
not vanishes). 

In the kinematical region under consideration, the divergence of the current (c) (see 
Eq. flFlSP ) is obtained by using Eq. flG6P : 

Q . {k'rk[Af''^krk..Yi, = ^ e eiirgf 6 ~ ^ " 



(K+ - k[^) A^i^y 



x^+ik'ic 



nil 
Pt 



itKiPlon) - 1 



F-A+(A;i,„)A+(A;^,jF2„A+(A;: 



2on 



(G16) 



The divergence of the current (i), Eq. (lF24p . in this region, using Eq. (1G8|) results 



A 



xA+(A;U)r? 



1 - A+ Pio„, — 7i 
Pi 



A+ ( A;io„) A+ (/C2o„)F2aA+ (/C2on) 



(G17) 



where we have used that Q = A^^^-* — A*^-*^-^^ + p'^^^ — k 
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Ion- 



Finally, adding Eqs. (IGlSp . (1G16P and (lG17p . only the second term in the square bracket 
of flG17|) survives, resulting in: 



X A+ ( A;;,jr? A+ ip[JA+ik'^jT2aA+ik2on)A+ (Ko„) ^7i+ A+ ( A^lon) = 

Pi 

= (A:i+^iJ^(i)(ir^)b;Vix)eiA+(p'ioJ^7i+A+(A:i,„,)^ {Q^ - k[^) . (G18) 

Pi 

This single term corresponds to the interaction calculated at first order in QP expansion, for 
total momentum Kf, in the considered kinematical region. The term with w^^\Ki) vanishes 
in this kinematical region, since we have Q'^ — k'^ = —p^ > and therefore we are outside 
the kinematical support of the interaction, see Eq. (IGlSp . 

The total divergence of the current operator is found by adding Eq. flG12p . flGlSp and 
Q . jc{o) ^ Then, by taking into account the definition of the left and right LP charge operator, 
Eqs. and (1661) . one can obtain an expression for the LP WTI valid for the whole 

kinematical range < A;^ < K'j . This illustrates how one can check the formal expression 
for Q-j'^^^\ given by Eq. (!74|) . in the example of the Yukawa model in ladder approximation. 
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